Microcavities coupled to multilevel atoms 
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A three-level atom in the A-configuration coupled to a microcavity is studied. The two transitions 
of the atom are assumed couple to different counterpropagating mode pairs in the cavity. We analyze 
the dynamics both, in the strong-coupling and the bad cavity limit. We find that compared to a two- 
level setup, the third atomic state and the additional control field modes crucially modify the system 
dynamics and enable more advanced control schemes. All results are explained using appropriate 
dressed state and eigenmode representations. As potential applications, we discuss optical switching 
and turnstile operations and detection of particles close to the resonator surface. 

PACS numbers: 42.50.Ar, 42.50.Pq, 42.60.Da 



I. INTRODUCTION 

Whispering gallery microcavities are an attractive im- 
plementation for quantum optical setups, often moti- 
vated by the very low loss rate that can be achieved. 
This has led to a multitude of applications p], , such 
as switches 0-0] , transistors H, 131, quantum networks 
[ToL 11 L photonic entanglement pjj], or optomechan- 
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The light is confined to the resonator by total inter- 
nal reflection, which leads to an exponentially damped 
evanescent field around the cavity. On the one hand this 
allows to couple the cavity, e.g., to an external fiber. But 
on the other hand, the evanescent field also allows the res- 
onator field to interact with nearby particles. Since the 
cavity volume and the resonator linewidth can be small, 
even the near-resonant coupling of a single atom can have 
a substantial effect on the resonator dynamics [HI and 
modify the photon statistics of transmitted photons . 
In [2(| , a photon turnstile was created by coupling a two 
level atom to the cavity. The sensitivity of the resonator 
to nearby atoms can also be used to measure the pres- 
ence of nearby particles, their optical properties, or their 
concentration |2l| - |27| . The evanescent field can also be 
used to trap atoms, as has been demonstrated for atoms 
coupled to a tapered fiber [28-3(1 ■ 

In most cases, the particle interacting with the res- 
onator is treated as a two-level atom, but there is some 
related previous work on more complex level schemes. 
For example, in Q a three level atom in A configuration 
coupled to a nearby microcavity was considered. Here, 
however, the atom only couples to the cavity with one of 
its transitions, while the other one is coupled to an ex- 
ternal laser field in order to control the atom-resonator 
interaction. In f3lj . a three level V system inside a c avity 
was considered to create entangled photon pairs. In [32|, 
a three-level atom in V-configuration coupled to a res- 
onator is studied, focusing on the effect of the scattering 



between the different cavity modes on the optical spec- 
tra. Transmission and reflection amplitudes for single 
photon wave packets interacting with a resonator coupled 
to a three- level atom are considered in [33|. Finally, in 
a related work [34j, photon scattering from a three-level 
atom coupled to a ID wave guide was studied. These 
works are not least motivated by the fact that a setup 
in which several transitions of the particle are coupled to 
respective modes of a single resonator, potentially could 
be more flexible and could enable more advanced control 
schemes. 

To demonstrate this, here, we analyze advanced control 
schemes in the generic model system of an atom coupled 
to a resonator, probed by a fiber coupled to the resonator. 
In particular, we study a three-level atom in the A config- 
uration coupled to a microcavity, see Fig.[T] We consider 
the case in which each atomic transition couples to a pair 
of counterpropagating cavity modes, respectively. In ad- 
dition, the resonator is coupled to a fiber which allows to 
probe the atom-cavity setup. We analyze the dynamics 
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Figure 1. (Color online) Schematic setup. The upper panel 
* sandra.schmid@mpi-hd.mpg.de shows the resonator coupled to the fiber and an atom. The 

t joerg.evers@mpi-hd.mpg.de lower panel shows the internal structure of the A-type atom. 
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Figure 2. (Color online) Alternative implementation of our 
model system based on two microcavities. Without the atom, 
the resonators do not interact, as they are not in resonance. 
The atom in A-configuration acts as a quantum link, with 
the two transitions each of which coupling to one of the res- 
onators, respectively. 



and the optical properties of the system in dependence on 
the fields applied via the fiber using two different approx- 
imations with complementary validity ranges. First, we 
assume weak input fields, and restrict the Hilbert space 
to low excitations. Second, we adiabatically eliminate 
the cavity and derive an effective master equation for the 
atomic degrees of freedom alone. We find that the exis- 
tence of two coupling mode pairs and two atomic ground 
states crucially changes the systems dynamics as com- 
pared to the two level case and opens up new possibilities 
to control the coupled system. For this, we typically treat 
one of the fields as probe field, and the other as control 
field. As potential applications, we discuss controllable 
optical switching and photon turnstile operations, as well 
as the measurement of the position of particles close to 
the resonator. All results are interpreted based on the 
system's dressed states or eigenmodes. 

Our system could be realized, for example, by tuning 
the Zeeman splitting of two atomic magnetic sublevels to 
the frequency spacing of two mode pairs in the resonator. 
Alternatively, different polarizations couple be used to 
selectively couple the transitions to different modes, but 
the leakage pattern of evanescent fields with different po- 
larizations around the resonator surface is not uniform. 
Finally, also a setup with two nearby resonators as shown 
in Fig. [2] could be considered in which the resonators do 
not interact in the absence of the atom as they are not in 
resonance with each other. The atom could then act as 
a quantum link between the resonators, with each tran- 
sition coupling to one of the resonators, respectively. 

This article is organized as follows. In Sec. Ill Al we 
introduce our model system derive the equations of mo- 
tion. In Sec. Ill B[ we describe the two calculation meth- 



ods based on a truncation of the Hilbert space (jll B 21) or 
on an adiabatic elimination of the resonator (jll B 3|) used 
throughout the analysis in detail. Our observables, the 
transmission and reflection intensities as well as second 
order correlations are introduced in Sec. IIIII Section IIVI 
discusses our results. We present our findings for the 
strong coupling regime in Sec. I IV Al and explain them 
via a dressed state atomic level scheme. Afterwards we 
show how our setup in this parameter regime could be 
used as a precision measurement device for the position 
of a nearby atom. In Sec. IIVBI we show our results for 
the bad cavity regime. We present transmission and re- 
flection intensities as well as our studies on the photon 
statistics. Finally in Sec. IIV B 31 we show how our cou- 
pled system can be used as a tunable bimodal photon 
turnstile in the bad cavity regime. 



II. THEORETICAL CONSIDERATIONS 

A. Model system 

The system we investigate is shown in Fig. [T] It con- 
sists of a toroidal microresonator, which is driven by a 
fiber, and coupled to an atom placed within the evanes- 
cent field surrounding the resonator. The atom is mod- 
eled with three states in the A-configuration with a en- 
ergy separation S of the two ground states. The res- 
onator is modeled with two pairs of counterpropagat- 
ing whispering gallery modes, with one pair coupling to 
each transition in the atom, respectively. This selective 
coupling could be achieved either via different mode fre- 
quencies, or via different polarizations. An alternative 
possible realization is shown in Fig. [3J and will be dis- 
cussed in more detail in Sec. |V] In our system a num- 
ber of different coupling mechanisms are of importance. 
First, within the resonator, each pair of counterpropa- 
gating modes with photon destruction operators a^, bi is 
coupled via the scattering rate hi (i € {1, 2}). In princi- 
ple, also photons from one mode pair could be scattered 
into the other pair of modes. This coupling is described 
by the parameters p and q. But in particular if the two 
mode pairs have different frequencies or polarizations, p 
and q will be negligible. Next, the fiber couples via the 
evanescent field to the cavity such that photons of the two 
probe fields with coupling strength E\ and £2 can enter 
the cavity. This coupling is described by the constants 
Krm.ex, where m € {a, b} and i G {1,2}. The incoherent 
decay of the cavity modes into the vacuum we include by 
the internal decay rate K mi ,in- Then, the overall loss rate 
for the cavity modes is 

^rrii ^mi,in ~t~ ^mi,ex • (1) 

Finally, the coupling of the atom to the resonator is de- 
scribed by the constants g mi . The respective coupling 
strength depends on the position of the atom character- 
ized by the radial position r and the azimuthal position 
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x along the resonator circumference, and is of the form 

9a.=9ti{r,z)e ik >\ (2a) 
9 bi = 9* at ■ (2b) 

Here, ki is the wave number of the mode pair i and z is the 
vertical coordinate. The first mode pair % — 1 couples to 
the first atomic transition |1) -H> |3), whereas the second 
pair i — 2 couples to the second transition |2) -f-> |3). We 
further define the atomic transition operators 



s? = m\ 
sr = m\ 



(3a) 
(3b) 



Since in general the frequencies of the driving fields, the 
resonator modes, and the atomic transition frequencies 
may be different, we introduce the detunings 



A; = (w 3 - u)i) - col , 



(4a) 
(4b) 
(4c) 



where uj 1 c are the two cavity mode pair resonance fre- 
quencies, uip the laser frequencies of the probe fields in 
the fiber, and W3 — Ui — cu l a the atomic transition fre- 
quencies. 

With these definitions, the Hamiltonian "Hat, for our 
system can be written as 

T~iab = Hab + ^ab + ^ab ; 



(5a) 

2 2 

i—l i—l m£a,b 

U s ab = hp(a\a 2 + b\b 2 )e iet + hq{a\b 2 + b\a 2 )e iet 
2 

+ h i a lh : + H.c. , (5c) 



n% = nY J £> l + hY J E 9* m M s i + Hc - . ( 5d ) 

i—l ra^a.b 



i=l 



where T-L° ab contains the free energies of the atom and the 
cavity, 1-L s ab the scattering between the different modes 
within the resonator, and l-i c ab the coupling between fiber 
and resonator as well as from resonator to the atom. The 
subindex ab indicates that the Hamiltonian is written in 
terms of the resonator modes di and bi. Including the 
incoherent processes leads to the master equation 



dtp = - r [T-Lab, p] + C K p + £ 7 P , 

2 

C K p =E E K m . (2rnipm\ - rn\mip 

i=l me{a,b} 



-pm\rrii 



(6) 



(7) 



2 

LrP = E I ( 25 r^+ - 5+5r p - pS+Sr) , (8) 

z=l 



where C K p describes cavity decay, and £ 7 p the atomic 
spontaneous emission. 

We next introduce the normal modes 13511 



B t = 



V2 ' 
ctj - bj 

V2 ' 



(9a) 
(9b) 



which leads to a partial diagonalization of the Hamilto- 
nian Eqs. ([5]) to give 



W-AB — 1~L AB + H C AB + H AB + 1L AB , 
2 2 

h°ab = - ft E A * 5 r 5 < + + ^E^ + 



i=l 
2 



2 



K AB = ^=Y { {S*A + £*B l )+ H.c., 



(10a) 

(10b) 
(10c) 



«is = ^E ^ 5 i" - H». B / 5 r) + H.c. , (lOd) 



H AB = h(p-q)e iet B+B 2 



h{q + p)e iet AfA 2 + H.c. 



(lOe) 



where = cos(fcix) and = sin(fcia;). 

Transforming also the incoherent parts to the new ba- 
sis, we obtain the equations of motion for the normal 
modes as 



d t Ai(t) = - T[Mt),U AB ] - K Ai At 
h 

£. 

= - i(Si + hi)Ai - i-j= - igA,S[ 
- i{q + p)e ltt A^ - K Ai Ai , 
d t B i {t) = - % -[B i {t),U AB ]-K Bi B i 



(11a) 



i(Si - hi)Bi - i— I - gBiSi 
i(p - q)e iet B^i - n Bi B t , 



(lib) 



where the two mode pairs are labeled by i € {1, 2}, and 
-^i is 1 for i = 2 and vice versa. Since modes of one 
pair {cii,bi} are of the same frequency, it is reasonable 
to assume K ai = K bi = Ki holds. The incoherent part 
remains diagonal also after the basis transformation, and 
we can write K Ai — K Bi — Ki and analogously for the 
internal and external decay rates K Auex — n Bi _ ex — Ki^ ex 
and KAi, in — K Bi,in = ^i.in, respectively. As expected, 
it can be seen from Eqs. (JXTJ) that for p = = g, the 
different normal modes Ai and Bi become independent 
apart from the coupling via the atom. 
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B. Two calculation methods 

In this section, we outline the approaches used in the 
following to solve the given system. We start by charac- 
terizing the two parameter regimes studied throughout 
the later analysis. Then, we describe two different meth- 
ods to obtain the steady state solution for our model sys- 
tem. In the first method, we assume weak input helds. 
This allows us to approximate the infinite series of Fock 
states of each photon mode by the lowest few photon 
number states. This truncation of the Hilbert space ren- 
ders the system finite, and thus allows us to numeri- 
cally calculate a steady state density matrix. We call 
this method the truncated Hilbert space (TH) method, 
and discuss it in Sec. Ill B 21 In the second model, we as- 
sume the bad cavity regime for the resonator, such that 
the loss and photon scattering channels of the resonator 
dominate the system dynamics. Then, the cavity modes 
can be adiabatically eliminated. This method will be 
termed adiabatic elimination (AE) method, and is ana- 
lyzed in Sec. lIIB3l Finally, in Sec. lIIB4l we compare the 
two methods and discuss their respective validity ranges. 

1. Parameter regimes of interest 

We consider two fundamentally different parameter 
regimes [36|. On the one hand, we study the strong 
coupling case, in which the coupling between the atom 
and the cavity dominates the dynamics compared to the 
loss dynamics of the cavity modes. Thus, the conditions 
Ki <C g mi are fulfilled. For the resonator-fiber part of the 
system, we assume critical coupling characterized by 

Ki,ex = ^? + «? 1<B . (12) 

As the rates Ki contain the coherent scattering rates hi 
via Ki >ex as well as the incoherent photon decay rates 
Ki t i n , in the strong coupling regime, we also have hi < 
Hi <C g mi ■ On the other hand, we consider the bad cav- 
ity regime. Here, the incoherent dynamics of the cav- 
ity modes dominates the dynamics, and the condition 
n m,i <C Ki holds. As before, then also g mi <C hi is ful- 
filled. 



2. TH method - Fock mode truncation 

The general state of our system is determined by the 
state of the atom, as well as the state of the four different 
cavity modes. The dynamics can thus be described using 
basis states of the form 

\atom, n(Ai), n(B 1 ),n(A 2 ), n(B 2 )> , (13) 

where n(A4) denotes the number of photons in mode M.. 
If the parameter are chosen such that the mean occu- 
pation number of a cavity mode remains low, then it is 



possible to restrict the corresponding state space to low 
Fock number states. Taking into account states with at 
most / photons in each mode would lead to a state space 
reduced to dimension 3 x (I + l) 4 , since also the state 
without photons has to be considered for each mode. In 
the following, we assume the case in which at most one 
photon can be in each of the cavity modes. The number 
of basis states then reduces to 3 x 2 4 = 48, and the cor- 
responding reduced density matrix pFock has 48 2 = 2304 
elements. Note that taking into account up to two ex- 
citations per mode would already lead to 59049 density 
matrix elements, which is impractically large even though 
the number can be considerably reduced due to symme- 
tries of the density matrix. 

We arrange all elements of pFock in a vector and 
eliminate the last element using the trace condition 
Tr(pFock) = 1- The resulting vector pFock has 48 2 — 1 
elements, and we can rewrite the equations of motion for 
the density matrix elements as 

dtPFock = Q ■ PFock + K . (14) 

The steady state pFock,stst can then be obtained by solv- 
ing the system of linear equations to give 

PFock, stst = —Q~ l ■ K . (15) 



3. AE method - Adiabatic elimination of the cavity modes 

In this approach, we assume that the dissipative dy- 
namics of the cavity dominates the system dynamics, 
such that the degrees of freedom of the cavity can be 
adiabatically eliminated from the system 37]. In partic- 
ular, we assume the bad cavity condition g mi <C Hi to 
be fulfilled. After the elimination, we obtain an effective 
master equation for the atomic degrees of freedom only. 
This significantly reduces the system dimension and thus 
simplifies the solution. 

We apply standard techniques to eliminate the res- 
onator modes [37[. Details on the calculation can be 
found in Appendix [SJ After the basis transformation 
Eqs. ©, we can treat the two normal mode pairs Ai and 
Bi separately. For each mode pair, first, we assume low 
occupation of the resonator modes due to the large dissi- 
pative dynamics, and expand the total density matrix in 
the photon occupation number. In the resulting density 
matrix equations, we assume 



Ki 

where |£ tom| represents the magnitude of terms corre- 
sponding to the coherent dynamics of the system, with 
| • | a suitable norm. Setting furthermore the time deriva- 
tives of off-diagonal elements in the expanded density 
matrix to zero, and inserting the thereby obtained ex- 
pressions for the off-diagonal elements into the diagonal 
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elements, we finally obtain an effective master equation 
for the atomic dynamics alone given by 



the strong coupling regime, only the TH method can be 
applied. 



OtQatom — ~ Jj[Hatom, £?atom] + £r£*atom , (17a) 

-pS+Sf) , (17b) 

2 2 
— 1 2=1 

+ MkSi + + Ml*Sr . (17c) 

The constants Ay and fi^ are effective detunings and 
Rabi frequencies arising from the cavity elimination, re- 
spectively. Their values depend on the various coupling 
constants governing our system's dynamics. Explicit ex- 
pressions can be found in Appendix [Bj Note that related 
results were obtained for a two-level atom coupled to a 
resonator with only one mode pair in [20| . 

4- Comparison of the two calculation methods 

The adiabatic elimination method leads to an effec- 
tive master equation, which can easily be solved and 
thus considerably reduces the computational effort. This 
comes at the price that the method only applies to the 
bad cavity limit, in which the dissipative dynamics dom- 
inates the system. In contrast, the Fock state truncation 
method does not a priori exclude either the bad cavity 
or the strong coupling regime, but rather puts a con- 
straint on the mean occupation number of the cavity 
modes alone. While potential resonant enhancement of 
the cavity occupation for low cavity loss rates must be 
taken into account, a weak input field generally suffices to 
remain within the validity range, such that this method 
also allows to access the strong coupling regime. Draw- 
backs of this method are the higher computational effort 
compared to the AE method, and the fact that by con- 
struction higher order correlation functions involving two 
photons are not accessible from the results. We further 
analyzed this point for a two-level atom coupled to the 
resonator, for which higher Fock modes can be included 
while keeping the total number of equations reasonably 
low. We found that it is not sufficient to take Fock states 
with two photonic excitations in total into account in 
order to reliably calculate the second order correlation 
function. Rather, one has to include two excitations per 
mode. As already mentioned, in the three-level atom 
case, including up to two photons per mode would in- 
crease the state space substantially, effectively rendering 
the calculations unfeasible. 

Thus in summary, in the bad cavity regime, both meth- 
ods can be used, and we verified that the two methods 
give the same results. Due to the reduced computational 
effort, here, the AE method will be used. In contrast, in 



C. Input-Output Relations 

In this section, we relate the fields inside the resonator 
to the fields in the fiber. First, we derive an expression 
for the driving terms proportional to Si in the equations 
of motion for the modes Ijlip. For this, we take the ex- 
pectation value of Eqs. (jlip . and obtain 



\d t Ai) = (d t A l ) S y S 7= 



71 



(18) 



and equivalently for modes Bi. The index sys indicates 
all internal terms of the atom-cavity system not related 
to the external driving. Note that we neglect the fluctu- 
ations of the input field in the driving terms entering the 
equations of motion for the modes by equating (Si) = £j. 
On the other hand, written in terms of an external in- 
put field operator Ai^ n , the corresponding equation for 
( d t Ai ) can be expressed as 

(d t Ai) = (d t A l ) S ys + y/K^- (Ai, m ) . (19) 

From Eq. (fl8|) and (fT9)) . we can thus conclude 

i£i 



(Ai^in) 



(20) 



Note that if no input field is applied to the modes 
then (A,in) = -^T"- 

Next, we consider the output from the resonator to the 
fiber. Using the input-output formalism HHHl], we find 



i,out — Q>i,in(fy \J ^^i^ex^i (^) 
= — (o-i.in) ~ a i,in{t) 

i,out — — bi.in(t) ~f~ y2's|ex^i(^) 

= - b' itin (t) + y/K — [Ai{t) ~ Bi(t)] 



(21a) 



(21b) 



where a' iin (t) and &■ in (t) are fluctuations, such that for 
example a^m = (ai,i n ) + aj irl (i). As long as we con- 
sider normal ordered output fields we can substitute the 
output operators according to [2(| 

ai,out = &io + aa<Sj~ + ot i2 S 2 ~ , (22a) 
bi,out = Pro + PiiS{- + fc 2 S 2 - . (22b) 

Here, the complex coefficients aio and an and a l2 are 
calculated as follows. We solve the equations of motion 
Eq. (jllbjl for the normal cavity modes Ai and Bi in the 
steady state. For the modes Ai the terms without an 
atomic operator we call o^o > the coefficient of S{~ in the 
equation for Ai we call an and the terms in front of 
operator S 2 ~ are a^. Analogously for the modes Bi we 
define the coefficients as /3,-q, flu and fia. 
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III. OBSERVABLES 
A. Transmission and reflection spectra 

The transmission spectra can be calculated from the 
normalized first order correlation functions as 14011 



7>„. = 



( m i,out m i,out) ( m iout m i,out} 



■ t 

, a i. out a i, out 



(23) 



for modes characterized by m £ {a, 6} and i 6 {1, 2}. 
The output photon operators can be related via Eqs. (|22|) 
to the solution of equations describing the cavity-atom 
system. l~F a . can be interpreted as the transmission in 
dependence of the input field frequencies, whereas Tp b . 
is the corresponding reflection. The transmission and re- 
flection can be calculated by both the AE and the TH 
methods introduced in Section lll B[ and are thus accessi- 
ble in both the strong coupling and the bad cavity regime. 



B. Photon statistics 

In order to investigate the photon statistics of the out- 
put photons we have to consider the normalized second 
order correlation functions at equal times. They can be 
calculated as kfj 



(2) 



( m { OM t TO Lui TO j, outrrii, out ) 
! m l out m i, out ) ( mt out m^ out ; 

,( 2 ) 



(24) 



If i = j, the function g F gives the photon statis- 
tics of a single mode. If the correlation function is one, 
then the photons obey Poissonian statistics. Higher or 
lower values indicate super- or sub-Poissonian statistics, 

(2) 

respectively. In case i j, the function g F is a 

cross correlation function for photons arriving at the de- 
tectors at the same time but in different modes rrii and 
rrij. As it was found in [2(| for a two- level atom cou- 
pled to a resonator, and as we also find in our anal- 
ysis, the obtained photon statistics in the critical cou- 
pling regime is predominantly determined by the atom 
because the resonator contribution is suppressed by de- 
structive interference in the critical coupling case. This 
suggests that sub-Poissonian statistics can be associated 
with anti-bunching, whereas super-Poissonian statistics 
can be connected to bunching. The second order corre- 

(2) 

lation functions g F can only have positive values 

(2) 

> 0. The slope of g F ' at time delay r — de- 

termines whether the photons exhibit bunching or anti- 
bunching. Bunching is indicated by a negative slope of 

(2) 

g F whereas anti-bunching occurs for positive slopes. 



For g F 



(2) 



= the slope can only be non-negative and 



therefore in case of g F 
have anti-bunching. 



(2) 



we can conclude that we 



Since for these second order correlation functions two 
photon processes are crucial, we can calculate them only 
by using the AE method introduced in Sec. Ill B 31 but 
not via the TH method including only up to one-photon 
Fock states. Thus we only present results in the photon 
statistics in the bad cavity regime. 



IV. RESULTS 

In this Section, we present the results obtained nu- 
merically for our observables. As the system has many 
degrees of freedom, we have to make a number of param- 
eter choices. Throughout the analysis, we assume that 
the coupling between different pairs of modes is weak 
and set q = p = 0. Furthermore, we consider the nat- 
ural case that the two modes belonging to one WGM 
pair have similar properties, which leads to a number 
of parameter choices. First, we set k^i — kb 1 = k>i 



and k^ 2 



K B2 



K2 where Ki — 



and 



Ki,ex = \J K-i i n + tif . Similarly, for our calculations we as- 
sume similar properties of the two pairs of modes. We as- 
sume the prefactors in the atom-cavity coupling strengths 
Eqs. (|2|) to be equal, i.e. <?q = = go. Note, how- 
ever, that the different position-dependence still in gen- 
eral leads to different couplings of mode Ai and Bi to 
the two atomic transitions. For the coupling between 
resonator and fiber we set K ex ,i = n e x,2 = n ex , for the 
atomic decay rate 71 = 72 = 7 and for the scattering 
inside the cavity hi = = h. The internal decay rate of 
the cavity modes are set to Hi^ n = 7, from which follows 
Ki = «2 = K. Furthermore, if not noted otherwise, we 
choose kiX = tt/2 which results in g^ = and g# i — go- 
Also, we assume that both the cavity mode pairs are on 
resonance with the atomic transition they couple to, re- 
spectively. Then, for the detuning Si — Ai holds. The 
energy separation between the two ground states is set 
to 8 = 0. 



A. Strong coupling regime 

In this section, we show our results for the transmis- 
sion and reflection in the strong coupling regime. For the 
calculation we use the TH method described in Sec. IIIBI 
As parameters we choose go = IOO7 and h = I57. First, 
we present our results for the output flux of the differ- 
ent modes. Afterwards, we explain them using a dressed 
state model of the coupled atom. 



1. Transmission and reflection 

Fig. [3] shows the transmission T~F ai and reflection 7~F bl 
in dependence on the two probe field detunings Ai and 
A2 In the following, we analyze different regions in this 
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Figure 3. (Color online) In (a) the normalized output flux 
7F b of mode bi tOU t in dependence of Ai and A2 is shown, 
(b) shows a cut of l~F ai and T~F bl at A2 = 1497, which is close 
to one of the eigenvalues of our system. The parameters are 
(50, h, Ki,in)h= (100,15,1). 



figure and interpret the results with the help of the eigen- 
states of the system and the populations in the different 
atomic states. For this, the relevant state space deter- 
mined below is shown in Fig. 01 and the populations of 
the atomic states in Fig. [5] 

If the absolute value of one of the detunings |Aj| is 
large, then the corresponding atomic transition \i) — > 
1 3) essentially decouples from the cavity. This leads to 
optical pumping into the atomic state such that the 
corresponding population approaches one. For example, 
at Ai = — 25O7 and IA2I < |Ai| the atom is more likely 
excited from state |2) to state |3) than it is from |1) to 
|3). Therefore, in the longtime limit the atom moves into 
state |1). Eventually, the atom can be neglected, and the 
two mode pairs decouple. As for large detuning the light 
passes the resonator unperturbed, the reflection Tf 6i is 
negligible. 

A special situation arises if Ai = A2. Then, both 
transitions of the atom are driven in the same way, and 
the atom evolves into a dark state [4l| with both ground 
states equally populated, and no population in the upper 
level V3, as shown in Fig. [5] This suppresses the interac- 
tion of the atom with the resonator, and thus the results 
on the main diagonal line Ai = A2 of Fig. [3] are equal to 
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Figure 4. (Color online) Dressed state atomic level splitting 
for excitation with up to one cavity photon per mode. 



the case without atom. 

The further structures of the atomic populations Vi 
and the reflection T~F b can be understood from an eigen- 
state analysis of the matrix governing the equations of 
motion of our system. For this, we diagonalize the sys- 
tem's internal Hamiltonian without external input field 
and thus calculate the eigenenergies and the dressed 
eigenstates of the coupled cavity-atom system. We note 
that our parameter choice KiX — 7r/2 leads to = and 
thus only the modes Bi can couple to the atom. Then, 
the relevant basis states reduce to \atom, n(Bi), 71(62)). 
Including up to one photon per cavity mode, five differ- 
ent eigenvalues and eleven eigenstates are obtained. The 
eigenvalues for the parameters used in Figs. [3] and [5] are 



Egl — Eg2 — . 
1 



E e2 = -h=-15>y 

E e4 = X ~ (-3/1 



8<4 + hp 1 
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(25a) 

-149 7 , (25b) 
(25c) 

134 7 , (25d) 



77 7 , 



(25e) 



E, 



e5 



-1227, (25f) 



and the corresponding eigenstates read 

\V g .i) = 1 1,0,0) , |^, 2 > = 12,0,0), 



(26a) 



\V el ) =|1,1,0) + |2,0,1) 



lE, 



CO 



9B l 



3,0,0), (26b) 



|y e2) i) = |2,l,0), |K 2 , 2 ) = |1, 0,1), (26c) 
|F e2>3 ) =-|l,l,0) + |2,0,l) 

\V e3 ) = |1,1,0> + |2,0,1>- — 13,0,0), (26d) 

= i(E e2 + 2h) + 

I V e ,, 2 ) = l{E&2+2h) % 1,1) + |3,1,0), (26f) 
9b z 
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Figure 5. (Color online) In (a) the population of the atomic 
state 1 is shown in dependence of both the detunings. The 
parameters are (go,h, Kj,m)/7 = (100, 15, 1). In (b) the pop- 
ulation of state 1 3) is shown. 
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We can identify three different kinds of eigenstates. The 
first group contains the ground states \V 3t x) and IV^), 
with the atom in one of its two ground states and no pho- 
ton in the cavity. Since we typically apply weak probe 
fields, our coupled system is mostly in one of these two 
states. The second group contains all states containing a 
single excitation of the system. The eigenstates belong- 
ing to this group are \V e \), | V^2,i— 3) and \V e z). These 
three states \V e 2) are the states in which the atom is in 
one of the ground states but the cavity is populated by 
one photon. The states \V e \) and \V e ^) are superposi- 
tions including also the state with excited atom and no 
photon in the cavity modes. The third group contains 
states in which the system contains two excitations, one 
of mode B\ and one of mode 82- It should be noted that 
these states cover only a small part of the doubly excited 
eigenstate space due to our initial choice of considering at 
most one photon per mode. But due to the low excitation 
considered, the doubly excited states have only negligible 
weak influence on the system dynamics and therefore can 
be neglected for the interpretation. The states including 



up to one excitation are shown in Fig. [4] 

To interpret the system, it is important to note that 
due to our definition of the detunings, the probe fields are 
resonant to an atomic dressed state transition if the de- 
tuning has the opposite sign of the respective eigenenergy 
shift. First, we focus on the line defined by A2 = 1497, 
for which the probe field £2 is resonant with the tran- 
sition \Vg^) — > \V e i). Fig. \3[b) shows the transmis- 
sion l~F ai and the reflection l~F bl at this detuning. For 
Ai = — 2507, the condition |A 2 | -C |Ai| holds and thus 
l~F bl ~ and the atom is almost completely in state |1), 
see Fig. [5j Scanning Ai towards zero detuning, a res- 
onance structure with increased reflection can be seen 
at around Ai = — 1347. At this point, the field £\ is 
near-resonant with the transition \V g .i) — > |V e 3), which 
leads to a scattering of photons in the reflection chan- 
nel via the atom. This can also be seen by the increase 
in the excited state population of the atom, see Fig. [5] 
Further increasing Ai, the reflection becomes minimal 
at Ai pa — 8O7. In this region, field E\ is relatively far 
detuned from the eigenstates of the system, as can be 
seen from the high population of the atomic state |1). 
As V3 is relatively small, also the reflection remains low. 
At Ai = — 157 = —h, a strong resonance is observe in 
the reflection spectrum. This peak can be attributed 
to the non-coupling cavity mode A\. A closer analy- 
sis shows that the structure around Ai = — 157 = — h 
in fact contains two separate peaks, which leads to the 
apparent asymmetry. The second peak is located at 
Ai = 157 = h = —E e 2 and is due to the coupling to 
the eigenstates \Ve.i). At positive detunings Ai, we first 
encounter a region with low reflection, which can be inter- 
preted in the same way as the region around Ai = — 8O7. 
At detunings around Ai = 1497, both probe fields cou- 
ple to the eigenstate on the transition \V g _\) — > \V e i). At 
exact two-photon resonance Ai = A2, a narrow dip in 
the reflection can be seen, which is due to the above- 
mentioned dark state. In this case the system behaves as 
if the atom was not present. In the vicinity of the dark 
state, the population of the excited state is relatively high 
due to the near-resonant coupling to the eigenstate \ V e i), 
which leads to the increase in reflection. 

Next, we analyze the region around Ai ss in 
Fig. EJa), where high reflectivities occur. It can be seen 
that this structure has a small residual dependence on 
A2. At detunings A2 for which the second field is in 
resonance with one of the system's eigenstates, the re- 
flection resonance is centered around Ai = — h = — 15 j. 
For all other values of A2, the structure is centered ap- 
proximately around Ai = 0. The reason for this is as fol- 
lows. If the second field is detuned from all resonances, 
then the atom is optically pumped into state \V g 2) by the 
first field, and the system again acts as if no atom was 
present. Without atom, the two cavity mode pairs are 
symmetrically split due to their coupling h such that they 
have resonances at ±h — ±157. As these two resonances 
are close to each other compared to the line width, they 
overlap and appear as a single resonance centered around 
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Ai = 0. In contrast, if A2 is tuned to resonance with one 
of the system's eigenstates, then the atom contributes to 
the system's dynamics. Then, in the reflection spectrum 
the non-coupling mode at A2 = — 157 dominates, as al- 
ready explained for the central resonance of Fig. [3jb). 
Considering the atomic population of the excites state Vz 
in Fig. [5] we observe a peak for Ai = at A2 = 1347 and 
A2 = — 1497 whereas for the third eigenvalue A2 = 167 
the population V3 ~ 0. This is because the eigenstates 
Ve2,i-3 do not include any state in which the atom is 
excited. 



2. Dependence on the atom position 

In this Section, we analyze the dependence of the re- 
sults on the position of the atom. In general, the position 
of the atom determines the coupling constants gj^ and 
gt3 i , which in turn influence the spectral properties of the 
total system. A potential realization for the following 
analysis could be an atom trapped in a potential created 
by the evanescent field of the resonator. This would lead 
to a trapping in radial direction, but could still allow for 
movements of the particle along the azimuthal direction. 

In Fig. |51 we show our results for A2 = for different 
positions of the A atom. Here, the distance to the cav- 
ity is constant but the azimuthal position along the res- 
onator circumference and thus the phase kiX is varied. If 
kiX = 7r/2, the atom couples to modes Bi only as assumed 
for the previously shown results. In case of kix = 0, the 
atom couples to the A4 modes only. Intermediate values 
lead to coupling to both modes. Changing kiX from to 
7r/2 corresponds to moving the atom over a distance of 
A/4 along the circumference. It can be seen in Fig. [B] that 
while the general structure of the spectrum remains con- 
stant, the positions of both sidebands oscillate around 
their mean position. Also the height of the sideband 
peaks changes, with maxima for the two cases in which 
a single mode pair couples to the atom. At kiX = it/A, 
the sideband peaks vanish. Additionally, the slight asym- 
metry in the central resonance around Ai = changes. 
For example, in moving from kiX — to kiX — 7r/2, the 
structure becomes mirrored at the Ai = axis. This is 
because the modes Ai and Bi exchange their roles as cou- 
pled and uncoupled modes. For kix — OAn, both mode 
pairs Ai and Bi couple to the atom, but with different 
coupling strengths. For kiX = 7r/4 the atom couples with 
equal strength to the Ai and Bi modes. Therefore, the 
middle peak structure is symmetric around Ai = 0. 

The radial dependence is directly governed by the ab- 
solute value of the coupling constants gj^ and gs t . The 
peaks in the spectrum lie at the eigenenergies given in 
Eq. (1231) . For fixed A2 they move to higher values of | Ai | 
when increasing the strength of the coupling between the 
atom and the cavity. 

This dependence of the spectral properties on the par- 
ticle position can also be used to detect the presence of 
the particle, or to measure its position. 
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Figure 6. (Color online) Reflection Tp b for different azimutal 
positions x of the atom and fixed A2 = 0. We chose kx — 0, 
f , 0.4tt, and f . 

B. Bad cavity regime 

In this Section we analyze the bad cavity regime, and 
choose go — 7O7 and h = 25O7, and all other parame- 
ters are chosen as in Sec. IIV AI Our observables are the 
output fluxes and the second order correlation functions. 
For the calculations we use the AE method as described 
in Sec. HTB31 

1. Transmission and reflection 

We start by analyzing the transmission into mode 
a>i,out shown in Fig. [7] As in the strong coupling case, 
the transmission approaches one if the first probe field 
is far detuned, i.e., for large |Ai|. Residual interactions 
lead to a trapping of the atom in state |1) if the sec- 
ond field is less detuned. Decreasing the absolute value 
of Ai, the transmission reduces and approaches zero al- 
most independently from the second detuning A2. There 
is a single resonance structure around Ai = — 197 and 
A2 = — 197 with increased transmission 7~F a ■ A cut 
through Fig. [Ffa) at A 2 = -227 is shown in Fig. Efb). 
A closer analysis revealed that both the position and the 
amplitude of this peak depend on the atom position via 
the phase kiX. In Fig. kiX = n/2 which means that 
modes Bi couple to the atom whereas modes Ai do not 
couple. On the other hand, choosing kiX = such that 
Ai couples and Bi not, the resonance structure moves 
to Ai = 197, A2 = 197, and has the same intensity as 
for kiX = 7r/2. For values of kiX between and ir/2 
the resonance moves along the diagonal line and has re- 
duced amplitude. For equal coupling of the two modes, 
9 A = SB) the resonance vanishes. We can thus directly 
trace the resonance back to the coupled atom. A similar 
resonance was also found for a two level atom coupled 
to a resonator in the strong coupling regime [20| . Inter- 
estingly, however, the resonance is not located at one of 



10 



(a) 



3001 



200 \ 



100 \ 



< j 

-iooj 

-200 1 



-300^ 



1O.6 



0.5 



0.4 



0.3 



0.2 



0.1 



(b) 



-400 -300 -200 -100 100 200 300 400 





1 1 1 







-200 -100 O 100 200 

Aj/7 

Figure 7. (Color online) (a) TF ai for the bad cavity regime 
with go = 707 and h = 25O7. (b) Cut at A 2 = -227 for the 
same parameters. 



the eigenstates of the system obtained by diagonalizing 
the coherent part of the initial interaction Hamiltonian. 
Rather, the resonance coincides with maximum popula- 
tion in the atom's upper state. This can be explained by 
diagonalizing the effective Hamiltonian obtained after the 
adiabatic elimination in Eq. I17cl The only nonvanishing 
eigenvalue of the system after the adiabatic elimination 
is equal to the sum of the corrected effective detunings 

An + A 22 . For our parameters An = A 22 = pr+^r- In- 
serting the parameters yields an expected position of the 
resonance of Ai = A 2 = —19. 57. Finally on the diago- 
nal line Ai = A 2 , the atom is in a dark state such that it 
decouples from the resonator. Thus the transmission is 
as in the case without atom for these parameter values. 
In the lower subfigure, the dark state leads to the sharp 
dip around Ai = —227. 



2. Second order correlation functions 

In this Section, we analyze the second order correla- 
tion functions g F . First, we consider the case i = j. 

In Fig. |8l we show the second order correlation func- 
tion for mode a\. ou t- In the lower subfigure, a cut at 
A 2 = —227 is shown. The dominant feature is a peak 
around Ai = I3O7 and A 2 = — 3O7. Around this posi- 



(2) 

tion, g F reaches high values up to 10. Thus, for these 
detunings the photons obey a super-Poissonian statistics. 
This structure is due to small values of Tpy for these 

(2) 

detunings. Since the correlation function g F is nor- 
malized to the output flux, see Eq. [Ml small intensities 

(2) 

lead to an enhancement of g F . Another important 
feature is the large area around zero detunings in which 
the system exhibits sub-Poissonian statistics. This area 
coincides with the resonance in the transmission around 
Ai = — 197, A 2 = — 197 found in Fig. [7J In this param- 
eter regime, non-classical light with sub-Poissonian pho- 
ton statistics is generated. Since this second order corre- 
lation function has values around zero, we can conclude 
from this information that we have anti-bunching. The 
non-classical light arises from the interaction with the 
atom, similar to the turnstile operation reported in [20| . 
Without atom, the transmission is negligible due to de- 
structive interference in forward direction downstream 
of the resonator. But excess excitations emitted from 
the atom can be transmitted in forward directions, and 
these exhibit sub-Poissonian statistics and anti-bunching 
due to the fact that the atom requires a finite time to 
be excited. Finally, in the dark state case Ai = A 2 , the 

(2) 

photon statistics is Poissonian, g F =1. 

Next, we considered the cross correlation functions. In 
Fig. O the cross correlation of mode ai >(m t an d a 2jOMt is 
presented. In can be seen that it also is close to zero in 
the region around Ai = —197, A 2 = — 197 which is the 
position of the transmission resonance in Tf q1 . This ex- 
tends the turnstile operation to two distinct modes ai. out 
and a 2 . ou t. Thus a detector placed in the transmission 
direction only detects individual photons in either of the 
two modes. We also observe two regions with high second 

(2) 

order correlation g F in Fig. [DJa). Analogously as for 

(2) 

g F , this can be explained by minima in either of the 
two output intensities to which the correlation function is 
normalized. In the region around (A 1, A 2 ) = (130,-30)7 
the first order correlation l~F ai ~ 0. In the second peak 
region around (Ai,A 2 ) = (—30,130)7 the transmission 
of mode a 2 almost vanishes. 



3. Applications 

In this Section, we outline possible applications of our 
setup in the bad cavity regime. First, we show that the 
coupled cavity atom system can be used as a single pho- 
ton source. Next, we show that it acts as a bi-modal pho- 
ton turnstile. Finally, we show that our bimodal photon 
turnstile is tunable by the input fields and their detun- 
ings. This could be exploited as an optical switch. 

Considering Fig. |SJa), we see that in a region of the 
resonance in TF ai close to (Ai,A 2 ) = (—19,-19)7 (c.f. 

Fig- 0) the second order correlation g F vanishes. This 
means that photons of mode a\. ou t leave the system indi- 
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Figure 8. (Color online) (a) Second order correlation function 
g F for the bad cavity regime with go = 70-y and h = 25O7. 
(b) Cut at A2 = —227 for the same parameters. 

vidually, such that it is impossible to detect two photons 
of this mode at the same time. Moreover, the cross cor- 
relation g F also vanishes, see Fig. 19(a). This means, 
the photons of each mode come out individually and ad- 
ditionally, if an di photons is detected no a% photon can 
be detected at the same time. Consequently, whenever an 
a>i,out photon is detected, it is impossible that an d2 l0 ut 
photons arrives at the detector at the same time. There- 
fore we know, that in a direction only individual photons 
leave the systems. The output intensity of both trans- 
mission modes on tOU t is different from zero Tf 01 7^ and 
7F a2 7^ 0, see Fig. [Jj Thus our coupled system can be 
used as bi-modal single photon source in a,i direction. 

This can also be used as a photon turnstile, in a di- 
rection. Photons come out individually and turn their 
propagation direction in case of non vanishing atom cav- 
ity interaction, which means that the atom is not in a 
kind of dark state. The turnstile can be controlled by 
changing the detunings and therewith driving the atom 
into a dark state. Controllable here means that we can 
turn the output flux on or off simply by changing the 
detuning of one of the input probe fields. For example, 
if we choose A2 = — 197 an d Ai w — 197 but 7^ — 197. 
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Figure 9. (Color online) (a) Cross correlation gp 1 for the 
bad cavity regime with go = 707 and h = 25O7. (b) Cut at 
A2 = —227 for the same parameters. 



Then the transmission intensities Tf and Tf are al- 
most maximum and thus 7^ 0. However, when slightly 
tuning the input laser field £\ to Ai — A2 = —197, 
the atom is in the dark state and thus does not interact 
with the cavity any more. Then the system behaves as 
if no atom was there and the transmission of both the 
modes is almost zero. Thus, whether we have a nonvan- 
ishing transmission intensity or not can be controlled by 
slightly changing either of the input detunings. 

Furthermore, an optical switch for the modes di jOU t 
can be realized by our system in the bad cavity regime. 
If the parameters are such that the atom is effectively 
decoupled from the resonator, due to the critical cou- 
pling condition no incoming aj,j n photon can be trans- 
mitted. Thus all incoming photons are reflected back 
from the cavity. But enabling the atom-cavity coupling, 
some photons are transmitted and thus the flux in a di- 
rection becomes non-zero. As an example, we found that 
for Ai rj -197 but Ai ^ -197 and A 2 = -227 the 
transmission Tp ai is maximum, see Fig. [7(b). But if the 
input field £2 is switched off, the atomic population is 
completely transferred to ground state |2) such that the 
atom decouples from the resonator. Then, the transmis- 
sion Tf„ ~ 0, and all light is reflected. 
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SUMMARY AND DISCUSSION 



and the atomic states 



We presented a system consisting of a tapered glass 
fiber and a toroidal microcavity coupling to a 3-level 
atom in A configuration. We considered two different pa- 
rameter regimes, the strong coupling regime and the bad 
cavity regime. For both these regimes we presented ade- 
quate calculation methods in order to solve the equations 
of motion for the system's operators, namely, one based 
on a truncation of the Hilbert space, and one based on 
an adiabatic elimination of the resonator. We compared 
the respective validity ranges, and their advantages and 
disadvantages. Our observables are the transmission and 
reflection output fluxes as well as second order correlation 
functions and cross correlations and the atomic popula- 
tions. We found that our system can work as a tunable 
bimodal photon turnstile which can be controlled by the 
detuning of the two input fields and explained how to 
exploit our setup as a photonic switch. Additionally, we 
analyze the dependence of our observables on the posi- 
tion of the nearby atom. In both the strong and weak 
coupling regime, we explained the structures of the trans- 
mission and reflection using respective dressed states or 
eigenmodes of the system. 

In our calculations, we neglected direct interactions 
between different mode pairs in the resonator, and as- 
sumed selective coupling of the two atomic transitions to 
one mode pair, respectively. This could be realized by 
using the similar frequency, but different polarizations 
for the two mode pairs which, however, is technically 
challenging. Alternatively, different frequencies and sim- 
ilar polarizations could be used, and the atom could be 
tuned in resonance with both modes by adjusting a Zee- 
man shift using an external magnetic field. Finally, our 
system can also be realized by two cavities coupling to 
the same A atom. This alternative implementation for 
the considered setup is shown in Fig. [2] In this setup, 
the atom couples to two resonators, and each mode pair 
is restricted to one of the resonators. As the modes in 
the two resonators have different polarization or different 
frequencies, direct coupling between the resonators could 
be avoided. This setup at the same time could provide a 
connection to studies on quantum transport of electrons 
through single molecules |42j |. 



ACKNOWLEDGMENTS 

Financial support by the International Max Planck Re- 
search School for Quantum Dynamics in Physics, Chem- 
istry and Biology is appreciated. 



Appendix A: Details on the adiabatic elimination 



d t p = -^[HAB,p] + ^ K p + C^p (Al) 
with the Liouville superoperators 
C kP = km(2M P M^ - M^Mp- pM^M) 

M£{Ai,B z } 

(A2) 

and 



= E ?( 25 r - S+S-p - pS+S-) . (A3) 



Due to the weak input field and the fast decay of the cav- 
ity population, we can assume that the average occupa- 
tion number is low. Furthermore, after the basis trans- 
formation Eq. ([9]), we can treat the two normal mode 
pairs Ai and Bi separately. We show the elimination 
formalism with modes Ai in the following. We start by 
expanding the density matrix in the photon number oc- 
cupation, keeping terms including up to one photon per 
mode. The expanded density operator reads |37j ] 
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(A4) 



where the photon states for the two modes of the Ai pair 
are defined as |n(^4.i), n(A?)). The pij are the residual 
part of the density matrix. 

Substituting the expanded density matrix Eq. (jA4|) 
into the master equation I All we obtain equations of mo- 
tion for the elements p^. The equations corresponding 
to the lowest photon number occupation evaluate to 

dfPoa =£atomP00 + ^\P\\ + 2^2^22 

1 



—f^2 + 9A 2 S 2 ) P20 



We start with the Hamiltonian given in Eq. (|10e|) and 
the complete master equation including all cavity modes 



- Pax +9A 1 S 1 

-P02 (~^^ 2 + 9A 2 S^ 
9tPW =£atomP00 — KlPlO 

- i( (5t + hi)p 10 



(A5a) 
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9A 2 S 2 + -j= s 2 ) P30 + {m+ p)e let p 20 ) , 



(A5b) 



dtPll —CatomPll — 2Kl/?ll 



1 



* ' ' TI^ 1 + 9AiSi i m 



-Pio +9AtS-L 
Here, C a tom.p is defined as 



(A5c) 



C-atomP — — -r[Hoi p] + £"iP i 
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2 



(A6) 
(A7) 



In the bad cavity regime, the cavity loss dynamics dom- 
inates over the coherent dynamics, such that 



(A8) 



«1, 



and thus the terms C a tom can be neglected in Eqs. (|A5|) . 
Setting further the time derivative of the diagonal ele- 
ments in Eqs. (|A5I) to zero, we arrive at 



Pw 



«i + i(Si + hi) 



+9AxSi)poO 

Pl2(-^£2 + 9A 2 S 2 ) 
'(9A 2 S 2 + + -^£ 2 )P30 + (m + p)e let p 20 



(A9) 



Now we calculate the time evolution of the atomic den- 
sity matrix using the diagonal elements of d t p exp 

dtgatom = 9 t poo + d t pu + d t p 22 . (AlO) 

Since our cavity is almost always empty, we substitute 
the terms obtained for the off diagonal elements pij only 
up to poo order. 



Therewith we finally arrive at an atomic master 
equation of the form as shown in Eq. (|17c|) . 

Appendix B: Constants for the Atomic Master 
Equation 

The constants used in order to calculate the atomic 
master equation Eq. (|17c[) are defined as follows: 

An =9aM X a)+9bM X b) 
A 22 =9 2 aMa) + 9 2 b Mb) 
Ai 2 =A 21 = 

fii =&Ai + ^Bi 

Q, 2 =£Ia 2 + 

with 

=Kn9Ai -j= + F*AP9At -|= 
1 l A 2 -Ka9a 2 — 7= + J~aP 9a 2 — 7= 

S2 Bl =iX B g Bl —= + iT B vg Bl —= 
£ 2 £\ 

^B 2 =^B9B 2 -j^ + i^B v *9B 2 ^J= 

Tii =2 ■(9 2 a ^(X A )+9bM X b))+ 7i 

T 22 =2-(9aMU) + 9bMb)) + 12 

Ti 2 =2 ■ (j9Ai9A2I&(Fa) + 9b 1 9b 2 ^{J 7 b)) 

r 2 i =r; 2 

p =(q+p)e let 

=(p-iy et 

X A = fM 

X B = 

U = 

& = 

Ta = 
Tb = 

fAi = 
fBi = 

!ai = 

fs 2 = 



1 - 


fAifA 2 PP* 




/Bi 


1 - 


IbJb 2 w* 




fA 2 


1 - 


!aiJa 2 pp* 




fs 2 


1 - 


!bJb 2 w* 




fAi Ia 2 


1 - 


!aJa 2 pp* 




/Bi fB 2 


1 - 


fB 1 fB 2 W* 

i 


KAi 


-i(5 1 + h 1 ) 
i 


«Bi 


- i(Si - hi) 
i 


KA 2 


- i(S 2 + h 2 ) 




i 



kb 2 - i{&2 - h 2 ) 



(BI) 
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